Eigenstates of a Small Josephson Junction Coupled to a Resonant Cavity 



W. A. Al-SaidiB and D. StrouJ 

Department 0} Physics, The Ohio State University, Columbus, Ohio 43210 

(February 1, 2008) 



We carry out a quantum-mechanical analysis of a small 
Josephson junction coupled to a single-mode resonant cavity. 
We find that the eigenstates of the combined junction-cavity 
system are strongly entangled only when the gate voltage ap- 
plied at one of the superconducting islands is tuned to cer- 
tain special values. One such value corresponds to the res- 
onant absorption of a single photon by Cooper pairs in the 
junction. Another special value corresponds to a two-photon 
absorption process. Near the single-photon resonant absorp- 
tion, the system is accurately described by a simplified model 
in which only the lowest two levels of the Josephson junction 
are retained in the Hamiltonian matrix. We noticed that this 
approximation does not work very well as the number of pho- 
tons in the resonator increases. Our system shows also the 
phenomenon of "collapse and revival" under suitable initial 
conditions, and our full numerical solution agrees with the 
two level approximation result. 

PACS numbers: 03.65. Ud, 03.67.Lx, 42.50.-p, 74.50.+r 



I. INTRODUCTION 

Circuits involving small Josephson junctions have the 
potential to behave as macroscopic two-level systems 
which can be externally controlled. Indeed, several re- 
cent experiments have demonstrated that such a system 
can be placed in a coherent superposition of two macro- 
scopic quantum states. The experiments have involved 
small superconducting loops, and also so-called Cooper 
pair boxes. In the former case, it was shown experimen- 
tally that the loop could exist in a coherent superpo- 
sition of two macroscopic states of different flux through 
the loop. In the case of the Cooper pair box, the sys- 
tem was placed into a superposition of two states having 
different numbers of Cooper pairs on one of the super- 
conducting islands f^-|(|. Such a coherent superposition 
had been proposed by several authors ^,|). 

One possible application of such two-state systems is as 
a quantum bit (qubit) for use in quantum computing 
At sufficiently low temperatures, this Josephson qubit 
will have little dissipation, and hence may be coherent for 
a reasonable length of time, as is required of a qubit. But 
in order for it to be potentially useful in computation, it 
must be possible to prepare entangled states of two such 
qubits - that is, states of the two qubits which cannot be 
expressed simply as product states of the two individual 
qubits. 

Such entangled states have now been created and ma- 
nipulated in systems involving atoms in a high-Q cav- 



ity pfj| |. Only recently has there been any theoreti- 
cal study of entangled states involving Josephson junc- 
tion circuits. Shnirman etal. [|ll| have considered low- 
capacitance Josephson junctions coupled to an LC-circuit 
in the two level approximation. Buisson and Hckking 
p^t have studied the states of a single Cooper pair box 
coupled to the photon mode in an electromagnetic res- 
onant cavity. Both and considered junctions 
such that the charging energy is much larger than the 
Josephson energy. Everitt etal. |l3| have carried out a 
similar investigation of the states of a small supercon- 
ducting quantum interference device (a SQUID - i. e., a 
small superconducting loop) coupled to a resonant cav- 
ity. In all of these studies, entangled states of the photon 
and the Josephson junction emerged naturally from the 
quantum-mechanical analysis. 

In the present paper, we carry out a full quantum- 
mechanical investigation of a model for a Josephson junc- 
tion coupled to a single-mode resonant photonic cav- 
ity. Our model Hamiltonian is, as we demonstrate here, 
equivalent to that studied in Ref. jL2|, via a canonical 
transformation. However, we consider a different param- 
eter region where the charging and Josephson energy of 
the junction are comparable in magnitude. Also we go 
beyond the analysis of that paper by including a large 
number of eigenstates of both the Josephson and the pho- 
ton Hamiltonian, in calculating the energy levels of the 
coupled system. When the system is at resonance, we 
find that the entangled states are indeed, as expected, 
well described by a basis of two states of the combined 
system. We also find a variety of entangled states, de- 
pending on the value of the control parameter (an offset 
voltage), including, at certain values, an eigenstate which 
is a superposition of two states in which the number of 
photons differs by two. We have also verified that our 
system shows the phenomenon of "collapse and revival" 
similar to that seen in quantum optics Our approach 
is similar to that applied by Everitt et al. |ll| to a SQUID 
ring in a resonant cavity, though our Hamiltonian is not 
equivalent. 

The remainder of this paper is organized as follows. In 
Section II, we give the model Hamiltonian, present sev- 
eral of its properties, and describe our method of solving 
it. The following section gives numerical results obtained 
from the model. The last section contains a concluding 
discussion, an estimate for some of the parameters used 
in our model, and a description of possible problems for 
future work. 
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II. MODEL 



A. Model Hamiltonian 



We consider an underdamped Josephson junction in a 
large-Q electromagnetic cavity which can support a single 
photon mode. The system is assumed to be described by 
the following model Hamiltonian: 



H — Wphc 



H.j + He- 



(1) 



tor potential takes the form A = \Jhc 2 juj V (a + a*) e, 
where e is the unit polarization vector of the cavity mode, 
and V is the volume of the cavity. Here we assumed that 
the junction dimensions are much smaller than the wave- 
length of the resonant cavity mode, so that the cavity 
electric field is approximately uniform within the junc- 
tion. Given this representation for A, the phase factor A 
can be written as A = (g/y/2)(a + aJ), where 
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(2) 



and £ii = e • £ is the effective length across the junction 
parallel to e. The definition of the Hamiltonian is com- 
pleted by writing down the commutation relations for 
the various operators. The only nonzero commutators 
are [a, a'] = 1, and [n, (ft] = —i. 

It is useful to make a cha nge of variables p — 
iy/hcu/2 (a) — a) ;q = ^Jfij (2 lo) (a> + a) . Then p and q 
satisfy the canonical commutation relation [p, q] = —ih, 
and 7i p hoton can be expressed in a form which makes 



clear that it represents the Hamiltonian of a harmonic 
oscillator: 



Hphoton = ^ (p 2 + u 2 q 2 ) 



(3) 



In terms of these new variables, the gauge-invariant phase 
difference is given by 7 = (ft — g y^uj/hq. 

The Hamiltonian H can now be expressed conveniently 
as the sum of three parts: 



and 



Here 7i p hoton is the Hamiltonian of the cavity mode, 
which we express in the form Hphoton = frw (a^a + 1/2) , where 
where and a are the usual photon creation and an- 
nihilation operators. Hj is the Josephson coupling en- 
ergy, which we write in the form Hj = — J cos 7, where 
J is the Josephson energy of the junction and 7 is the 
gauge-invariant phase across the junction (defined more 
precisely below). J is related to I c , the critical cur- 
rent of the junction, by J = M c /(2\e\). Finally, He 
is the capacitive energy of the junction, which we write 
as He = \V ' (n — n) 2 , where U — 4e 2 /C is the charging 
energy of the junction with capacitance C; n is an oper- 
ator which represents the difference between the number 
of Cooper pairs on the two superconducting islands which 
make up the junction; and n is a tunable experimental 
parameter related to the gate voltage which is applied at 
one of the superconducting islands. 

The gauge-invariant phase difference 7 is the term 
which couples the Josephson junction to the cavity. It 
may be written as 7 = 4>— (27r/$o) / A-d£ = 4>—A, where 
$0 = hc/(2e) is the flux quantum, (ft is the phase differ- 
ence across the junction in a particular gauge, A is the 
vector potential, and the line integral is taken across the 
junction. We assume that A arises from the electromag- 
netic field of the cavity normal mode. In Gaussian units, 
and assuming the Coulomb gauge (V • A = 0), this vec- 



H — HjJ + Wphoton + Hint, 



Hjj = - U (n ~ n) 2 - J cos < 



Hint = — J (COS 7 — COS ( 



(4) 



(5) 



(6) 



In the absence of H m t the Hamiltonian is the sum of 
two independent parts, corresponding to the Josephson 
junction and the resonant photon mode. H m t couples 
these two together. 



B. Method of Solution 

It is convenient to diagonalize H in a complete basis 
consisting of the direct product of the eigenfunctions of 
Hjj and Wphoton- 

The eigenfunctions of 7i p hoton are, of course, harmonic 
oscillator eigenstates; the n th eigenstate has wave func- 
tion 



h n (q) 



\/v / ^2"n! 



eM-y 2 /m n (y), 



(7) 



where H n (y) is a Hermite polynomial of order n, and y = 
(ui/hY^q. An eigenfunction ift((ft) of Hjj with eigenvalue 
Ej(n) can be written as ift((ft) = exp(i n (ft) rj((ft), and 
satisfies the Schrodinger equation, Hjjift((ft) = Ejift((ft). 
This equation can be written out explicitly using the rep- 
resentation n = —i- 
tion relation between n and 



■i-TT, which follows from the commuta- 



a.s 



d 2 Y(x) ( {8Ej 



2Qcos2a: Y(x) = 0, 



(8) 



where x = (ft/2, Y(x) = rj((ft/2), and Q = 4 J/U. This is a 
Mathieu equation with characteristic value a = 8Ej/U, 
and potential of strength Q. The eigenvalues Ej are 
determined by the requirement that the wave function, 
ift{(ft), should be single valued, i. e., that ip((ft) = ip((ft+2Tr), 
or equivalently Y(x + ir) = exp(— 2inir) Y(x). We there- 
fore define Yk(x) as a Floquet solution of the Mathieu 
equation (^), with Floquet exponent v = 2 k — 2 ft, where 
k = 0,±1,±2, . . .. We denote the corresponding eigen- 
value of Hjj by Ej t k(n)- For < n < 0.5, the lowest 
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eigenvalue corresponds to k = 0, followed in order by 
k = 1, —1, 2, —2, . . .. The eigenvalues of Hjj are periodic 
in n with period unity, and are also symmetric about 
fi = 1/2. 

Finally, any eigenstate ^(4>, q), of the Hamiltonian 7i 
can be expressed as a linear combination of product wave 
functions consisting of eigenfunctions of 7ijj and 7i p hoton: 



*(</>, 9) = ^2A kn ip k (<p) h n (q), 



(9) 



where ipk(4>) = exp(in0)Yfc(</>/2), and Afc„ are the ex- 
pansion coefficients. The only term in TL which is not 
diagonal in this product basis is the interaction term 
7ij nt . The eigenfunctions and eigenvalues of the Hamil- 
tonian of Eq. (Q) follow from the Schrodinger equation, 
HV{<t>,q) = E^(d>, q). 



C. Canonical Transformation 

Before proceeding to the solutions of this Schrodinger 
equation, we show that Tl is equivalent, via a canonical 
transformation, to another Hamiltonian which has some- 
times been discussed in the literature pl| , ^2| . A similar 
transformation has been demonstrated in Ref. [jllj for the 
case of one and two junctions. The required transforma- 
tion is as follows: 



cj>' = (f>- g^/uj/hq 

n' = n 

p' = p + gVKuj n 

q' = q. 



(10) 



It is readily found that the new operators satisfy the com- 
mutation relations [(f>',n'] — i, [q',p'] = ih, with all other 
commutators vanishing. The Hamiltonian 7i can readily 
be expressed in terms of the new variables. The result is 



H' = \u'[n' 



) - J cos <p + - (p +uj q ) 

2 „2 



gvhiv pn + 



U hum g 
2(U + hug 2 )' 



(11) 



where n' = U n/(U + hug 2 ) and U' = U + hu g 2 . 

If we interpret the photonic mode as an LC resonator, 
then this Hamiltonian describes a junction which is ca- 
pacitively coupled to that mode. Obviously, this trans- 
formed Hamiltonian will have the same spectrum of 
eigenvalues as the original one. Although it is not ob- 
vious from the form of Eq. (|Tl|), we have numerically 
confirmed that the spectrum is symmetric with respect 
to n — 1/2, as is implied by the original Hamiltonian. 



III. NUMERICAL RESULTS 

We have diagonalized Ti [Eq. (||)] numerically, using 
the basis discussed in Section IIB. We used a truncated 



basis which is composed of the direct product of the low- 
est five energy levels of the junction (k = 0, ±1,±2), 
and the lowest ten of the electromagnetic field (n = 
0, 1, . . . , 9). We have also confirmed numerically that in- 
creasing the number of states in this basis has little effect 
on the eigenvalues at least up to E = 6 U . 

Fig. 1 shows the calculated eigenvalues of Tt, plotted 
as a function of n up to an energy of U, for coupling con- 
stants g = 0, g — 0.15 and g = 1.5. In all three cases, we 
have arbitrarily chosen Q = 0.7 and hui/U = 0.3, corre- 
sponding to a case where the charging energy, Josephson 
energy, and quantum h cu of radiation energy are all com- 
parable in magnitude. The degeneracy of some of the 
energy levels in Fig. 1(a) is broken in Fig. 1(b) (the nu- 
merical value of the second and third energy levels split- 
ting at n = 0.258 is given in the caption); the degeneracy 
is more noticeably broken in Fig. 1(c), where the coupling 
constant g is much larger. The degeneracy-breaking is, 
of course, caused by the perturbation Hint- Moreover, 
the magnitude of the energy splitting also tends to in- 
crease with increasing g. This behavior is a character- 
istic feature of "level repulsion" predicted by standard 
degenerate perturbation theory in quantum mechanics. 

Next, we discuss the time evolution of the system de- 
scribed by Hamiltonian (^), given its state at time t = 0. 
The expectation value of any operator O at time t is 



(0(t)) = Tr(p(t)G), 



(12) 



where p{t) is the density matrix of the system at time 
t. p(t) is given by p(t) = U (t) p(0)W (t) , where U(t) = 
cxp(—iHt/h) is the evolution operator, and p(0) is de- 
termined by the initial state of the system. The time 
average of the operator O is then given by 



«0» 



lim 



Tr(p(t)0)dt, 



(13) 



where the inner and outer set of triangular brackets de- 
note, respectively, a quantum mechanical and a time av- 
erage. 

As an illustration, we have calculated ((7ijj)) and 
((^photon))- For each operator, we carried out the calcu- 
lation making the arbitrary assumption that the state 
of the system at time t = is \a = i\/2; k). Here 
\a = is a coherent state of the electromagnetic 

field, i. e., an eigenstate of the annihilation operator, a 
such that a | a) = a\a). (We have confirmed numerically 
that the time averages are unchanged if the initial state 
is an eigenstate of the number operator a^a.) \k) is an 
eigenstate of 7ijj with quantum number fc, and we have 
considered two different initial states, corresponding to 
k = and k = 1. The density matrix p(0) at time t = 
is easily calculated once the initial state is specified; in 
obvious notation it is p(0) = \a = iV2; k)(a = i\/2; k\. 
We have calculated the time evolution using the param- 
eters Q = 0.7, hui/U = 0.3, and g = 0.15. 

The time-averaged energy ((Hjj)) for these parameters 
is shown in Fig. 2a for k = and k — 1 . Each curve shows 
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a strong structure near n — 0.26 (and n — 0.74). At these 
values of n, according to Fig. 1, the difference in energy 
between the ground state energy, and the first excited 
energy is ~ 0.3J7 = h~u>. The structure thus corresponds, 
we believe, to a resonant absorption of a photon by a 
Cooper pair in the junction - that is, at this value of n, 
Cooper pairs of electrons move from the ground state, 
\k = 0), of the junction to the first excited state, |fc = 
1). The corresponding structure for k = 1 (dotted line) 
arises from stimulated emission, in which a Cooper pair 
falls from the excited state to the ground state of the 
junction with the emission of a photon. The exchange 
of energy between the junction and the resonator can 
also be seen in Fig. 2b, where we show the time-averaged 
energy contained in the photon mode, ((7i P hoton)), for 
the same two initial states. ((Wphoton)) for both states is 
close to the unperturbed value ((n) + l/2)hu> = 5 hw/2, 
which is the energy contained in the photon field for the 
coherent state \a — iy/2), but one increases sharply, while 
the other decreases, near the resonance at n = 0.26 (and 
n = 0.74). 

The curves in Fig. 2 also show weaker structure near 
n w 0.06 ( and near 0.94 ). This structure corresponds 
to a process in which two photons are absorbed by the 
junction Cooper pairs. Specifically, the junction is ex- 
cited from its ground state to its second excited state at 
this value of n, with a consequent loss of two photons 
from the electromagnetic field. This is discussed further 
below. 

We have also studied the time evolution of this system 
under conditions such that n is fixed at a value of 0.258, 
near the principal resonance mentioned above. To see the 
time evolution, we prepared the system at time t = in 
the state |n = 1 ; k = 0), for which the photon resonator 
is in state n = 1 and the junction in state k = 0. We 
then allowed the system to evolve in time according to 
the evolution operator U(t) = exp(—iHt/h). In Fig. 3 
we show the resulting time-dependent expectation value 
of (Ti-jj), and of (Hphoton)- Also shown in the same fig- 
ure are the time-dependent probabilities for finding the 
junction in state k = 1 at time t and the resonator in 
state n — 1 at time t, given this initial state. Clearly, the 
system is oscillating between the states \n = 1 ; k = 0) 
and \n = ; fc = 1). This energy exchange is periodic in 
time: both (Tijj) and (7i p hoton) vary periodically in time, 
with a period T\ w 675 h/U for our particular choice of 
Q/U and Tllo/U . In the first half of this period, Cooper 
pairs absorb energy from the electromagnetic field and 
are driven into the first excited state, while in the sec- 
ond half, they lose energy to the electromagnetic field and 
fall back into the ground state. In the present model, this 
process continues indefinitely. In a more realistic model 
which includes dissipation, these oscillations would grad- 
ually decay in time. 

Note also that the existence of a periodic oscillation 
between one state and another is not dependent on the 
choice of the initial state. For example, if the initial state 
is chosen as \n — ; k = 1), the junction would interact 



with the vacuum fluctuations of the electromagnetic field, 
and would periodically visit the state \n= 1 ; k — 0). We 
have confirmed this behavior numerically in our model. 

In Fig. 4, we show (Hjj) and (7i p hoton) as a function 
of time for two different initial states: \n = 1 ; k = 0) 
(full line) and |n = 1 ; k = 1) (dashed line). (The first 
of these is also shown in Fig. 3.) Evidently both states 
oscillate periodically in time, but with slightly different 
periods. 

Fig. 5 shows the analogous time evolution of a system 
prepared in a fixed initial state at n = 0.061, where we 
expect to see a signature of a two-photon absorption. 
In this case, we prepared the system in the initial state 
\n = 2; k = 0). The system oscillates between the k = 
(ground state) and k = — 1 (second excited state) of the 
Josephson junction, which is separated from the ground 
state by an energy 2 h lo at this value of n. The number of 
photons in the system also oscillates periodically in time, 
between the values n — 2 and n = 0, with period T2 ~ 
7500 h/U. This period is acd bout eleven times longer 
than the period at n = 0.258 (single-photon absorption), 
because the matrix element connecting the two states at 
n = 0.06 is much smaller. 



IV. DISCUSSION 

The problem studied here closely resembles that of a 
two-level atom interacting quantum-mechanically with a 
monochromatic electromagnetic field. Specifically, when 
n is such that the field is approximately in resonance with 
a transition between the lowest and first excited state of 
the Josephson junction, then it is a good approximation 
to consider only four states in the basis. We now show 
that this approximation gives results in good agreement 
with those obtained using the full direct product basis. 

We consider a basis constructed out of the two lowest 
energy states of the junction, i. e., for k = , 1 and of 
the electromagnetic field, for n — 0,1. Using a four 
dimensional basis formed from the direct product of these 
states, \n = 0, 1 ; k = 0,1), we can express the total 
Hamiltonian of Eq. (^) as a four times four dimensional 
matrix. If we retain in the basis only the two states 
1 ; 1) and |1 ; 0) which are approximately degenerate, the 
Hamiltonian can be further reduced to the 2x2 matrix 

W = eo,i|0; 1>( ; 1| + e 1)0 |1 ; 0)(1 ; 0| 

+ 6 |0; 1)(1; 0|+£r|l; 0)(0; 1|. (14) 

The diagonal matrix elements take the values eo,i = 
E JA (0.258) + hoj/2; e h0 = £J J>0 (0.258) + 3 hu/2, where 
Ej t k(n) is the k th eigenenergy of Hjj at n. £1 is given 
by £1 = (0, l|iW int |l,0), where H int is given by Eq. 
(|6|). To first order in the coupling strength g, 7ii n t ~ 
— Jg\J (oj/K) q sin <fi, and hence, to the same order in g, 
£l w —(Jg/y/8) (fc = 0| sin0|fc = 1). A similar two-level 
approximation can be carried out for coupling between 
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the states |n; 1) and \n — 1; 0); in this case, the cou- 
pling energy is found to be £i yfn. However, as expected, 
we have found that the two-level approximation becomes 
progressively more inaccurate as n is increased. 

To illustrate this approximation, we consider the case 
where the energy levels eo,i and e^o are degenerate. Diag- 
onalizing the reduced Hamiltonian matrix, Eq. (|l4|), gives 
us the two energy levels E± — eo ± |£i|. These levels cor- 
respond to two entangled photon- Josephson states. The 
Rabi period of oscillation between them is T\ — nTi/\£i\. 
Using the parameters and wave functions corresponding 
to our calculation of Fig. 3, we obtain T\ as 667 Ti/U, 
which is in excellent agreement with the value found from 
the full numerical solution (675 fi/U), based on the com- 
plete basis and shown in Fig. 3. We conclude that indeed 
the reduced Hamiltonian, defined in Eq. (|l4|), is an ex- 
cellent description of the system, at least when the two 
lowest levels are close to resonance. 

The two-photon process discussed in Fig. 6 can also be 
understood using the reduced Hamiltonian: 

H = e 2 . \2; )(2; 0|+e ,-i|0; -1)(0; -1| 

+ 6 |2;0)<0;-1| + £||0;-1)<2;0|. (15) 

Here we consider only the two degenerate levels with en- 
ergy, e 2 ,o = £7 J|0 (0.061)+5»w/2; e ,_i = (0.061) + 

hu/2. The coupling strength £2 is given by £2 = 
(2,0|±W int |0,-l), where U mt is given by Eq. §). 
To second order in the coupling strength g, 7ii nt as 
— (1/2) Jg 2 (to/fi) q 2 cos(/>, and hence, to the same order 
in g, £ 2 ~ ~(Jg 2 a/2/8) (k = -l|cos0|fc = 0). Using the 
same parameters as before, we find that the Rabi oscil- 
lation period is T 2 — 7r7i / 1 ~ 5800 fi/U, compared to 
the value 7500 h/U obtained from the exact numerical 
solution. Thus the two-level approximation is not quite 
as good as in the one-photon case, at least for these pa- 
rameters. 

As mentioned in the Introduction, our Hamiltonian 
model shows the phenomenon of "collapse and revival," 
provided that the resonator is prepared initially in a co- 
herent state |lj]. This behavior arises from the inter- 
ference between the different Rabi oscillations which ex- 
ist when the initial state is a superposition of different 
photon number states. In Fig. 6. we show the time de- 
pendent probability that the junction is in the excited 
state, k = 1, when the initial state of the system is 
\a = iy/3;k = 0). Also on the same figure we show 
(dotted line ) the same probability calculated using the 
reduced Hamiltonian of Eq. (@ Q: 

v(t) = ---e-^ 2 j2 m r cos(2 ^ v ~ st) - (16) 

The agreement between the two results is very good, es- 
pecially for short times. For longer times or for larger a, 
the agreement between the two level approximation and 
the full numerical solution is less good, primarily because 
the number of photons in the cavity becomes large. 



Before concluding, we make some order-of-magnitude 
estimates for some of the parameters entering our model, 
and compare them with those in plausible junctions and 
cavities. First, the parameter g may be estimated using 
Eq. (J2J) . If we arbitrarily assume a cylindrical cavity of 
radius r and length d, the frequency of the lowest mode 
is approximately 2c/r, and may be either a TE or TM 
mode, depending on the ratio d/r. Substituting into Eq. 
(§), we obtain g ~ t\\s/8/ (137 rd) ~ 0.24t\\/Vrd, where 
e 2 /(hc) as 1/137 is the fine structure constant. To esti- 
mate huj/U, we use the same expression for lu as 2c/r 
and U = 4e 2 /C, where C is the junction capacitance. 
Taking C = x, where x is a characteristic length, we 
obtain Tilu/U ~ 137 x/r. Thus, a value of 0.3 for this 
parameter implies x/r ~ 0.002 or less. For a capacitance 
of 10~ 16 F, readily achievable within current technology 
Q, we have x ~ 10~ 3 cm, or U ~ 1 X 10~ 15 erg - 8 K, 
and hence r ~ 0.5 cm. The associated frequency would 
be lu ~ 100 GHz, which is in the range used in some 
recent experiments |H| . Of course the quantum of radia- 
tion is limited by the superconducting gap A of the super- 
conductor used. We also note that our choice Q = 0.7 
corresponds to J ~ 3 K for the above parameters, or 
equivalently, a critical current I c of about 0.1 /zA. The 
value of J should be of the same order as a quantum of 
radiation, hu>. 

Using the above estimates, and assuming a realistic 
junction with an insulating layer of thickness £ ~ 100 nm, 
the coupling strength g will be g ~ 10~ 5 if r ~ d ~ 
0.5 cm. If we make <i<r, as in a disk-like microcavity, 
we can presumably increase g substantially, but it will 
still remain much below the value of 0.15 used in our 
calculations. Such a value of g = 0.15, or even larger, is 
preferable because this would make the Rabi period of the 
same order or less than the coherence time in Josephson 
junctions. In a high-Q cavity, the dissipation-induced de- 
coherence time may be as large as 2 ns || , much smaller 
than would be produced byg~10 -5 . A larger g would 
decrease the Rabi time, and technological improvements 
will no doubt increase the junction coherence time, mak- 
ing these times more comparable. Another way to in- 
crease the effective coupling, and hence decrease the Rabi 
period, would be to use many identical junctions located 
in a region small compared to a wavelength. The effec- 
tive coupling is then found to increase as the square root 
of the number of junctions in the cavity On the ba- 
sis of all these estimates, we conclude that the two-level 
approximation tested in this paper should be excellent 
for the range of parameters likely to be achievable exper- 
imentally. Of course, the precise relevant parameters for 
a given experiment depend on the details of the geome- 
try, which have only been very crudely estimated in the 
above discussion. 

It is of interest to compare the present quantum prob- 
lem to another system which bears some similarities, 
namely an array of Josephson junctions placed in a res- 
onant cavity |15| . This system involves large numbers of 
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junctions, each coupled to a single-mode cavity such as 
is considered here. While the appropriate Hamiltonian 
for treating such an array resembles that of the present 
problem, there are important differences. First, the ob- 
served array behavior involves large numbers of photons 
in the microcavity, whereas the present model is suitable 
to the case of one or a few photons. Secondly, the IV 
characteristics of the array involve junctions which are 
biased onto either the resistive or the superconducting 
branch, whereas the present model considers only the lat- 
ter branch. Finally, on the basis of the above estimates, 
the junctions to which the present calculation applies are 
probably substantially smaller, with smaller critical cur- 
rents, than those studied in Ref. |lfj. Thus, a rather dif- 
ferent approach, such as those described in Refs. (Tt| Q 
is called for. 

To summarize, we have considered the interaction be- 
tween a small Josephson junction and the photon field 
in a resonant cavity, in the limit when dissipation can 
be neglected. We find that there is a strong interaction 
between the junction and the resonant photon mode at 
special values of the gate voltage, as defined by the vari- 
able n. At such values of n, the lowest two eigenstates 
of the system correspond to entangled states involving 
states of both the photon field and the Josephson junc- 
tion. The Rabi period of oscillation between these states 
is inversely proportional to a certain coupling strength 
g. We also found evidence for a two-photon absorption 
process at another resonant value of the variable n. Fi- 
nally, we found that the lowest entangled eigenstates are 
accurately described by a simple two-level model which 
is a truncation of the full model Hamiltonian, Eq. 

We conclude with a brief discussion of possible im- 
provements in the present calculations. First, of course, 
the effects of dissipation need to be included. In prin- 
ciple such dissipation can be included by coupling the 
Hamiltonian degrees of freedom to a bath of harmonic 
oscillators. If the density of states of this bath has the 
right frequency-dependence, this coupling leads to Ohmic 
damping p2|-p4[. Secondly, it would be of great interest 
to extend this work to more than one junction coupled 
to the same resonant cavity [|l6] . If the two- level approx- 
imation continues to be accurate in this case, the study 
of such a group of Josephson junctions would be greatly 
simplified. 

V. ACKNOWLEDGMENTS 

We are grateful for useful discussions with E. Almaas. 
This work has been supported by the National Science 
Foundation, through Grants DMR97-31511 and DMR01- 
04987. 

FIGURE CAPTIONS 



1. Lowest eigenvalues ElU of the system described 
by the Hamiltonian (|6j) , calculated as a function of 
n for the parameters hu/U = 0.3 and Q — 0.7, 
and (a) g = 0, (b) g = 0.15, and (c) g = 1.5. 
All eigenvalues are shown up to E/U = 1. Note 
in particular the breaking of the degeneracy of the 
second and third energy levels around n = 0.26, 
seen clearly in Fig. lc. The separation AE between 
these levels at n = 0.258 is AE/U = 0, 0.01, and 
0.06 in (a), (b) and (c) respectively. 

2. Time and quantum-mechanical average of the op- 
erator Hjj [Fig. 2(a)], and of H pho ton [Fig. 2(b)], 
plotted as a function of n for two different assump- 
tions about the initial state of the system: \a — 
iy/2; k = 0) (solid curves); and \a — k = 1) 
(dashed curves). The states are defined explicitly 
in the text. In both cases, we used the parameter 
values hu/U = 0.3, Q = 0.7, and g = 0.15. 

3. (a) Time-dependent value of (7ijj) (full curve) and 
of (Hphoton) (dashed curve) for h = 0.258. The 
system is prepared at time t — in state \n = 
1; k = 0). The other parameters are Q = 0.7, 
hu/U = 0.3, and g = 0.15. (b) Same as (a) ex- 
cept that we show the time-dependent probability 
for finding the Josephson junction (full curve) and 
the photonic resonator (dashed curve) in their first 
excited states. 

4. (a) Time-dependent value of (7ijj) for two different 
initial states: \n = 1 ; k = 0) (full curve), and |n = 
1 ; k = 1) (dashed curve). In both cases, n = 0.258, 
Q = 0.7, hu/U = 0.3, and g = 0.15. (b). Same as 
(a) except we plot (7Y p hoton). 

5. Same as Fig. 3, except that fi is fixed at 0.061, 
and the system is prepared initially in the state 
| rt = 2 ; k = 0) at time t = 0. Other parameters 
are the same as in Fig. 3. 

6. Time dependent probability that the junction is in 
the first excited state given that initially the res- 
onator is in the coherent state \a = iy/3) and the 
junction in the ground state. Also shown (dotted 
line) is the probability V(t) [Eq. ([l6|)] calculated 
using the two level approximation. Here we used 
n = 0.258, Q = 0.7, and hu/U = 0.3. 
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